## Gilbert lower bound

from PyM import *

def lb_gilbert(n,d,q=2):
    return int(ceil(q**n/volume(n,d-1,q)))
    

show(lb_gilbert(8,3))
show(lb_gilbert(9,3))
show(lb_gilbert(10,3))
show(lb_gilbert(11,3))
show(lb_gilbert(11,3,3))
show(lb_gilbert(21,5))